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ABSTRACT 

We present theoretical predictions of UV continuum luminosity function (UV LF) and Lya equivalent width 
(EW) distribution of Lyman alpha emitters (LAEs) in the framework of the hierarchical clustering model of 
galaxy formation. The model parameters about LAEs were determined by fitting to the observed Lya LF at 
z = 5.7 in our previous study, and the fit indicates that extinction of Lya photons by dust is significantly less 
effective than that of UV continuum photons, implying clumpy dust distribution in interstellar medium. We 
then compare the predictions about UV LFs and EW distributions with a variety of observations at z '--^3-6, 
allowing no more free parameters and paying careful attention to the selection conditions of LAEs in each 
survey. We find that the predicted UV LFs and EW distributions are in nice agreement with observed data, 
and especially, our model naturally reproduces the existence of large EW LAEs (> 240 A) without introducing 
Pop III stars or top-heavy initial mass function. We show that both the stellar population (young age and low 
metallicity) and extinction by clumpy dust are the keys to reproduce large EW LAEs. The evidence of EW 
enhancement by clumpy dust is further strengthened by the quantitative agreement between our model and 
recent observations about a positive correlation between EW and extinction. The observed trend that brighter 
LAEs in UV continuum tend to have smaller mean EW is also reproduced, and the clumpy dust is playing an 
important role again for this trend. We suggested in our previous study that the transmission of intergalactic 
medium for Lya emission rapidly decreases from z 6 to 7 by the fitting to Lya LFs, and this evidence is 
quantitatively strengthened by the comparison with the UV LF and EW distribution at z ~ 6.6. 
Subject headings: galaxies: evolution — galaxies: formation — galaxies: high-redshift — methods: numerical 



1. INTRODUCTION 

Detecting redshifted Lya emission with narrow-band imag- 
ing is a powerful strategy to seek for galaxies in the early uni- 
verse. Indeed, many galaxies have been detected through this 
method in the last decade (e.g., Cowie & Hu 1998; Rhoads 
et al. 2000; Taniguchi et al. 2005; Murayama et al. 2007; 
Ouchi et al. 2008), and they are called Lya emitters (LAEs). ' 
Because of their strong Lya emission lines, LAEs can be de- 
tected even at very high redshifts, as demonstrated by the fact 
that the spectroscopically confirmed highest-redshift galaxy 
so far (z = 6.96) was found by this method (lye et al. 2006; 
Ota et al. 2008). Moreover, they are an invaluable population 
to probe the cosmic reionization history because the strength 
and profile of Lya emissions from LAEs could significantly 
be affected because of absorption by neutral hydrogen in the 
intergalactic medium (IGM; e.g., Malhotra & Rhoads 2004; 
Santos 2004; Haiman & Cen 2005; Dijkstra et al. 2007a, b; 
Dayal et al. 2008; Mesinger & Furlanetto 2008). 

The physical properties (such as mass, age or metallicity) of 
LAEs and the connection with other high-z galaxy population 
[e.g., Lyman-break galaxies (LBGs)] have been poorly under- 
Electronic address: mark.kobayashi@nao.ac.jp 

' In this paper, so-called Lya blobs that extend on super-galactic scale are 
not considered, and they are treated as a separate population from LAEs that 
we focus on. 



Stood because of the faintness of their continua. On the other 
hand, their statistical properties such as luminosity functions 
(LFs) in terms of Lya emission (Lya LF) and rest-frame UV 
continuum luminosities (UV LF), and Lya equivalent width 
(EW) distributions, have more firmly been established be- 
cause of the increase of survey fields and available samples 
by different authors (Hu et al. 2004; Kashikawa et al. 2006; 
Shimasaku et al. 2006; Dawson et al. 2007; Gronwall et al. 
2007; Murayama et al. 2007; Ouchi et al. 2008). 

In our previous study (Kobayashi et al. 2007, hereafter 
KTN07), we have constructed a new theoretical model for the 
Lya LF of LAEs in the framework of hierarchical galaxy for- 
mation. It is based on one of the latest semi-analytic model 
for galaxy formation, the Mitaka model (Nagashima & Yoshii 
2004; see also Nagashima et al. 2005), in which galaxies 
are formed based on the standard structure formation the- 
ory driven by cold dark matter There are several theoretical 
models for Lya LF of LAEs based on analytic models (e.g., 
Mao et al. 2007; Dayal et al. 2008; Samui et al. 2009), 
semi-analytic galaxy formation models (e.g., Le Delliou et 
al. 2006; Orsi et al. 2008) or cosmological hydrodynamic 
simulations (e.g.. Barton et al. 2004; Nagamine et al. 2008; 
Dayal et al. 2009). The intrinsic production rate of Lya pho- 
tons within a galaxy is expected to be proportional to ionizing 
luminosity and hence it can be calculated if the star forma- 
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tion history (SFH) is known. Therefore the escape fraction of 
Lya photons from the galaxy is the key to predict LAE statis- 
tics. However, the escape fraction is difficuh to predict the- 
oretically for a realistic configuration of interstellar gas and 
dust, although there are some theoretical studies to predict it 
through radiative transfer models in a simplified geometry of 
interstellar medium (ISM) and dust (e.g., Hansen & Oh 2005; 
Verhamme et al. 2008). Hence, most models simply assume a 
constant escape fraction, which clearly contradicts the recent 
observational results from local star-forming galaxies (e.g., 
Atek et al. 2008; OstHn et al. 2009). The KTN07 model 
is unique because the two physical effects are incorporated 
in calculating the escape fraction: extinction by interstellar 
dust but with an amount that is different from that for UV 
continuum, and galaxy-scale outflows induced as supernova 
feedbacks. In KTN07, we have shown that our outflow+dust 
model reproduces the observational Lya LPs of the LAEs at 
z~3-6. 

UV LPs and distributions of Lya EW also provide impor- 
tant statistical information for LAEs in addition to Lya LPs. 
In particular, because Lya EW is highly sensitive to SFH, age, 
and metallicity as well as initial mass function (IMP), the ob- 
served EW distribution can provide an opportunity to test the 
validity of the LAE models. The purpose of this paper is to 
present a detailed and comprehensive comparison between the 
KTN07 model and the available observations of these statis- 
tical quantities at various redshifts. We will then focus on the 
two interesting issues described below from our analysis. 

The first issue is the existence of very large Lya EW LAEs. 
Theoretical models of stellar evolution predict that the max- 
imum Lya EW powered by star-formation activity with the 
Salpeter IMP and the solar metallicity is 240 A (e.g.. Chariot 
& Pall 1993; Schaerer 2003). However, observations have re- 
vealed that some fractions of LAE candidates have higher Lya 
EWs than the maximum (e.g., Malhotra & Rhoads 2002; Shi- 
masaku et al. 2006; Dawson et al. 2007; Gronwall et al. 2007; 
Ouchi et al. 2008). Although the existence of such high-EW 
LAEs has not yet firmly been confirmed because of faint UV 
continuum flux with large flux errors, such large EW LAEs 
could have significant implications for galaxy formation the- 
ory. In fact, these observational results have often been used 
to argue that some fraction of metal-free stars (so-called Pop 
III stars) and/or a top-heavy IMP are required (e.g., Dijkstra 
& Wyithe 2007). Another possibility is the enhancement of 
EW because of selective extinction for continuum photons by 
dust in dense clouds (Neufeld 1991; Hansen & Oh 2006), and 
recent observations give some supports to this interpretation 
(Pinkelstein et al. 2008, 2009a, 2009b). Since our model in- 
corporates the effect of extinction by dust, it is interesting to 
see whether such an effect can quantitatively explain the large 
EW LAEs in our model. 

The second issue is the trend of LAE EW distributions as 
a function of UV luminosity. In the rest-frame UV absolute 
magnitude Muv versus EWlyq plane, there is a lack of large 
EW LAEs with large UV luminosities, and the maximum EW 
systematically decreases with increasing UV luminosity (Shi- 
masaku et al. 2006; Stanway et al. 2007; Deharveng et al. 
2008; Ouchi et al. 2008). This trend in the M(UV)-EWLya 
plane was first reported for LBGs at z ^ 5 - 6 by Ando et al. 
(2006), and we here call this feature as the Ando effect. The 
physical origin of the Ando effect is not well understood. A 
similar trend has akeady been known for broad emission lines 
of active galactic nuclei, which is called as the Baldwin effect 



(Baldwin 1977), but the physical origin can be completely dif- 
ferent for the Ando effect on galaxies. Ando et al. (2006) sug- 
gested that this relation may be attributed to higher metallici- 
ties in the UV-bright galaxies. Another possible interpretation 
proposed by Ouchi et al. (2008) is that the average stellar pop- 
ulation of the UV-bright galaxies is older than that of the UV- 
faint galaxies. Schaerer & Verhamme (2008) and Verhamme 
et al. (2008) also provided some qualitative predictions to the 
origin of the Ando effect based on the results from their ra- 
diative transfer model, but it was not examined quantitatively 
whether their prediction reproduces observational distribution 
in M(UV)-EWLya plane. On the other hand, if the extinction 
effect by dust is important for the large EW LAEs as men- 
tioned above, extinction should also be relevant to the Ando 
effect. We will try to give a new theoretical explanation for 
the Ando effect based on our model. 

KTN07 found that, although the Lya LPs of LAEs can be 
reproduced by the model at z ^ 3-6, the model overproduces 
Lya LPs compared with observations at z > 6, because of the 
rather sudden decrease of the observed Lya LPs at z > 6. 
KTN07 suggested that this can be interpreted by a rapid in- 
crease of IGM opacity against Lya photons at z > 6, giving an 
interesting implication for the cosmic reionization. Here we 
revisit this issue in terms of the increased statistical quantities 
of UV LPs and EW distributions, and examine the KTN07's 
interpretation. 

The paper will be organized as follows: In §|2l we describe 
our theoretical model, especially for the UV luminosity and 
EW calculation. We compare the model results with the ob- 
served LAE statistical quantities at various redshifts 3 < z < 6 
in § [3] and discuss about the above two issues. After some 
discussions including implications for reionization (§111), the 
summary will be given in § |5] The background cosmology 
adopted in this paper is the standard ACDM model: JIm = 0.3, 
Qa = 0.7, Qb = 0.04, h = 0.7, and erg = 0.9. AH magnitudes are 
expressed in the AB system, and all Lya EW values in this 
paper are in the rest-frame. 

2. MODEL DESCRIPTION 

A detailed description of our outflow+dust model about 
Lya emission from galaxies in the Mitaka model is given in 
KTN07. Here we briefly summarize the essential treatments 
about the modeling of Lya emission, and some changes and 
updates from the KTN07 model. 

2.1. Lya Photon Production 

We consider only Lya photons produced by star formation 
activity. The Lya photon production rate in a star-forming 
galaxy is calculated from ionizing UV photon (A < 912 A) lu- 
minosity assuming that all the ionizing photons are absorbed 
by ionization of hydrogen within the galaxy, and Lya photons 
are produced by the case B recombination. The escape frac- 
tion of ionizing photons, flf^, from a galaxy is not exactly 
zero in reality, but it is generally believed to be small (e.g., 

Inoue et al. 2006). If ftZ'^ <^ 1, the assumption of ftZ'^ = in 
this work is reasonable because the Lya luminosity is propor- 
tional to (1 -ftsc'^). 

The ionizing photon luminosity is calculated using the stel- 
lar evolution model of Schaerer (2003, hereafter S03) assum- 
ing the Salpeter IMP in 0.1-100 Mq. Its metaflicity depen- 
dence is taken into account using the S03 model in a range of 
Z/Zq = 0-2. The gas metallicity of each model galaxy is cal- 
culated in the Mitaka model. In the quiescently star forming 
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galaxies, which do not experience a major merger of galax- 
ies, we relate the ionizing luminosity simply to star forma- 
tion rate (SFR), because they are forming stars approximately 
constantly and the mean stellar age weighted by ionizing lu- 
minosity at the time of observation is much smaller than the 
time scale of SFR evolution. However, in starburst galaxies 
triggered by major mergers of galaxies, SFR changes with a 
very short time scale and hence we exactly calculate ionizing 
photon luminosity by integrating the SFH. 

2.2. Lya Escape Fraction 
The observed Lya luminosity is determined by the es- 
cape fraction, ftZ", of produced Lya photons from their host 
galaxy. It is difficult to predict /i^c" of each model galaxy 
from the first principles in its realistic geometry of ISM and 
interstellar dust because of the resonant scattering of Lya pho- 
tons. The photon path could be very complicated depending 
on dumpiness and velocity fields of ISM, and extinction by 
interstellar dust for Lya photons could be different from that 
for UV continuum photons. We introduce a simple model 

where /i^c" is separated into two parts: 



.Lya r l-eXp(-T^^^°) 
/esc - /{) 



Lya 



(1) 



where fo represents the reduction of escaping photons by 
physical effects not caused by interstellar dust, and the rest 
of the r.h.s. is for extinction by dust in the slab geometry 
[Rybicki & Lightman 1979, eq. (1.30); see also Clemens & 
Alexander 2004, eqs.(l)-(3)]. 

There are some physical effects that may reduce /o from the 
unity. One such effect is too many times of resonant scatter- 
ings of Lya photons. When the scattering length is too small 
in some star forming regions in a galaxy, Lya photons will 
diffuse out only by random-walk process, and the escape time 
scale can be extremely long, effectively resulting in a small fo. 
It should be noted that a possible deviation of fl'Z'^ from the 
assumed value of ftl,^ = is also absorbed within fo. Another 
effect on fo is the absorption by neutral hydrogen in IGM. Lya 
photons that are blueshifted than the rest-frame of their host 
galaxy will significantly be attenuated by this effect at z > 3, 
while those redshifted will escape freely. We assume that this 
effect is independent of redshift at z < 6, because we expect 
that the evolution of absorption is not significant if the Lya 
line profile is similar for LAEs at different redshifts. How- 
ever, at the very high redshifts of z > 6 reaching the epoch of 
reionization, the damping wing effect of IGM absorption by 
the increase of IGM neutral fraction may become important, 
and in this case the damping wing would erase most of Lya 
photons including those redshifted than the rest-frame of the 
host galaxy, resulting in a rapid evolution of IGM absorption 
at z > 6 (e.g.. Fan et al. 2006). Therefore we treat the ef- 
fect of damping wing separately from fo, by introducing the 
IGM transmission T^y^, which is the fraction of Lya photons 
transmitted after the absorption effect of the damping wing, 
as done in KTN07. We emphasize that T^^^ represents only 
the effect of the damping wing at z > 6, and hence 7]^™ = 1 
at z ^ 6. The modest effect of absorption by IGM at z < 6 is 
effectively included in the parameter /q. The case of T^^^ < 1 



at z > 6 will be discussed in § 14.11 

The parameter t^^" is the effective optical depth of extinc- 
tion by interstellar dust for Lya photons. Here, the slab type 



geometry has been assumed for the dust effect rather than the 
screen type geometry. KTN07 found that the difference be- 
tween the models assuming the slab or screen geometries is 
negligible about LAE Lya LF predictions, and we here adopt 
the slab geometry because it has akeady been adopted for the 
dust extinction of continuum photons in the Mitaka model. 
In the Mitaka model, the dust optical depth for continuum 
photons are assumed to be proportional to the metal column 
density of cold gas, A'coidZcoia, and we also assume this propor- 
tionality for the dust opacity for Lya photons, t^^" . The metal 
column density can be calculated by the gas mass, galaxy size, 
and metallicity in the Mitaka model. Therefore, t^^" is given 
by 



Lya 



^ Dold^cold 
(A'coldZcold)o*^" 



(2) 



with a model parameter (A^coidZcoid)o^" that controls the 
strength of Lya photon extinction. As mentioned above, this 
parameter can be different from that for UV continuum pho- 
tons around 1216 A. We determine this parameter by fitting 
to the observed Lya LF of LAEs at z = 5.7 (Shimasaku et al. 
2006) independently of continuum extinction. 

We also consider galaxy-scale outflow by supernova feed- 
back as a potential effect that could enhance the escape frac- 
tion because the velocity difference of ambient ISM produced 
by outflow may greatly increase the Lya scattering length. 
The outflow is considered only for starburst population trig- 
gered by major mergers of galaxies, and outflow occurs by 
supernova feedback. We identify outflow-phase galaxies in a 
way consistent with the treatment of supernova feedback in 
the Mitaka model. Hence, there are four phases of galaxies in 
our model: quiescent (non-starburst), pre-outflow starbursts, 
outflow starbursts, and post-outflow starbursts. Although the 
distinction of model galaxies into these four phases seems 
somewhat arbitrary, we have done this based on physical con- 
siderations. Briefly, the outflow phase onsets when the energy 
injected to ISM by supernova feedback exceeds the binding 
energy of ISM gas in the galaxy halo, and it continues dur- 
ing the dynamical time scale of halo (^ r^/Vc = fesc, where 
Te and Vc are the effective radius of galaxy and the circular 
velocity of its host halo, respectively). The pre- and post- 
outflow starbursts are defined as galaxies before and after this 
outflow phase (see KTN07 for more details). The galactic 
wind terminates star formation, and then galaxies in the out- 
flow or post-outflow phases are passively evolving. We then 
introduce f™'^ instead of eq. ([T]) as ftic' in the model galax- 
ies under the outflow phase. This is motivated by a physical 

consideration that ftlc" becomes less sensitive to t^^^" in such 
outflowing condition because dust in ISM becomes sparse by 
the outflow and outflow drastically reduces the scattering op- 
tical depth of Lya by neutral hydrogen in ISM. Finally, we as- 
sume that galaxies in post-outflow phase do not produce any 
Lya emission, because ionizing luminosity should be reduced 
by the termination of star formation, and there is little amount 
of interstellar gas to absorb ionizing photons and produce Lya 
photons within the galaxies. 

2.3. LAE Model Parameter Determination 

Consequently, there are three model parameters about 

LAEs: (A'coidZcoid)o^", /o, and We assume that they 

are independent of redshift, because these parameters reflect 
the physics about Lya photon escape, that is independent 
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of redshift. These have been determined in KTN07 by the 
fit to the LAE Lya luminosity function at z = 5.7 measured 
by Shimasaku et al. (2006). We found a unique set of the 
best-fit parameters^, which are: (A^coidZcoid)!)^" = 8.0;!;} 4 x 
lO^i [Zq cm-2], fo = 0.23!j]:g^ and /q™"* = 0.36. The model 
parameters that are not related to LAEs in the Mitaka model 
are kept at the original values that have been determined by 
fits to observations of the local galaxies (Nagashima & Yoshii 
2004). As demonstrated by KTN07, the predictions by this 
model are in overall agreement with observations of LAE Lya 
LFs at various redshifts at z ^ 3-6. It may be rather surpris- 
ing, considering the complicated physics about Lya emission, 
that this simple phenomenological model can explain obser- 
vations with just three free parameters. The outflow-phase 
escape fraction f™^ is not much different from fo, and we 
get a reasonable agreement with observations even if we set 
ywind _ indicating that the outflow effect is not significant 
about the Lya escape fraction. This also means that the ef- 
fective number of free parameters is further reduced, i.e., just 
two. 

We will then compare the UV continuum luminosities and 
EWs predicted by this model with observations. It should be 
noted that all the model parameters have been determined by 
Nagashima & Yoshii (2004) and KTN07, and there is no free 
parameter that can be adjusted to the new data compared in 
this work. Therefore, the comparison with observations pre- 
sented below provides an objective test for the validity of the 
KTN07 framework of LAE modeling. 

2.4. Rest-Frame UV Luminosity 

For the purpose of this paper, we need to calculate UV con- 
tinuum luminosity of galaxies. We present UV LFs at the 
rest-frame wavelength of A =1500 A. The unabsorbed (intrin- 
sic) UV luminosity at this wavelength is calculated by using 
the S03 model in a similar way to the ionizing luminosity. 
The observable UV luminosity is then calculated taking into 
account dust extinction. The amount of extinction magnitude 
A\ for continuum photons as a function of rest-frame wave- 
length A has been calculated in the original Mitaka model, 
assuming that the optical depth tJ(A) is proportional to the 
metal column density of cold gas. The wavelength depen- 
dence of tJ(A) is determined by the Galactic extinction curve 
(Pei 1992), and Aa is calculated from tJ(A) assuming the slab 
type geometry, i.e., 

,r.-OAA,_ '^-^W[-r'dW] ... 

The proportionality constant between and metal column 
density has been determined in the Mitaka model to fit the 
observed local galaxies (Nagashima & Yoshii 2004). 

For the starburst galaxies, we expect that their dust opaci- 
ties are gradually reduced because SN explosions and subse- 
quent galactic wind would heat and remove their interstellar 
cold gas. It is difficult to predict analytically the amount of 
ISM left during the starburst activity. Here we simply assume 
that ISM decreases exponentially with an e-folding time of 
fesc around t = fwind- We tested the sensitivity of our result to 

^ These values are slightly different from those listed in Table 1 of KTN07, 
because the dust geometry has been changed into the slab model in the 
dust+outflow model, and there was a tiny bug in ionizing luminosity calcu- 
lation of the previous analysis. However, the overall agreement between the 
KTN07 LAE Lya LF model and the observations is not significantly affected. 



this prescription by using another model simply assuming no 
extinction in the outflow and post-outflow galaxies, and we 
confirmed that our main conclusions are not significantly af- 
fected. 

It is interesting to compare the extinction of UV continuum 
photons to that for Lya photons in our model. Both t^" and 
rJ(A) are assumed to be proportional to the metal column den- 
sity of ISM, and the relation between the two around the Lya 
wavelength becomes 

T^^" =qdT^{XLya), (4) 

where = 0. 149 ± 0.03 is derived from the model parameters 
in our model. We note here that is not a new free param- 
eter in our model, because both t^^" and Tj(\i,ya) have been 
determined by the modeling described above. The numerical 
constant in eq. ^ is called as the geometry parameter or 
dumpiness parameter introduced originally by Finkelstein et 
al. (2008).^ The value of qd effectively reflects the interstel- 
lar dust geometry. The case of ^ 1 means homogeneous 
ISM, in which the resonance scatterings make the photon path 
of Lya photons much longer than that of UV continuum pho- 
tons and hence Lya photons suffer from much larger extinc- 
tion. On the other hand, the case of q^ <C 1 is regarded as 
extremely clumpy ISM, in which Lya photons are reflected 
by neutral hydrogen before entering the dense regions where 
a large amount of dust particles reside. In such a case qd be- 
comes smaller than unity, since Lya photons propagate pref- 
erentially in low extinction regions, while U V continuum pho- 
tons go through the dense dusty regions. Detailed theoreti- 
cal studies solving radiative transfer of Lya photons in such 
clumpy ISM found that this is indeed possible (Neufeld 1 99 1 ; 
Hansen & Oh 2006). The value of = 0.15 derived from our 
best-fit model implies strong dumpiness of interstellar dust 
distributions in high-z LAEs. 

2.5. Rest-Frame Lya Equivalent Width 

Next we calculate Lya EW of each galaxy in the model de- 
scribed above. In addition to the observable EW (i.e., dust 
extinction incorporated into both luminosities of Lya line and 
UV continuum) denoted as EWlyq, it is convenient to define 
intrinsic EW"y^ that is simply calculated using the S03 stellar 
spectra and the Salpeter IMF, assuming the case B recombi- 
nation and 100% escape fractions both for Lya and UV con- 
tinuum photons. Note that all EW values quoted in this paper 
are in rest frame. We need to calculate UV continuum lumi- 
nosity at A = ALya = 1216 A to estimate EW. We calculate it 
from the UV luminosity at 1500 A and assuming a spectral in- 
dex of (3 = -2, where f\ cx A^ (i.e., a flat spectrum in /^), for 
the stellar spectra without extinction by dust. This is a typical 
UV spectrum of young stellar population (Bruzual & Chariot 
1993; Leitherer et al. 1999), and it is often assumed in studies 
of high-redshift galaxies including LAEs (e.g., Ouchi et al. 
2008). The time evolutions of EWj^^^ in the two cases of in- 
stantaneous starburst population and constant star formation 
are plotted in Fig.[T](left panel), for some values of metallici- 
ties. The maximum values of EW"y'^ are obtained at an age of 
^ 1 Myr after the onset of star formation, with the values of 
- 240, 420, and 820 A for metallicities of Z/Zg = 1, 1 /2000, 
and (Pop III), respectively. 

' However, our definition of is slightly different from the original ob- 
servational definition of ^^^s by Finkelstein et al. See § 12.51 
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However, because of the Lya photon escape fraction and 
extinction of continuum photons, EWlyq can be different 
from EWj"y^. From the modeling above, the ratio is given 
as: 

^ _ EWLyg ^ /O l-eXp[-g^T^(ALya)] 

~EWi";„"^rf l-exp[-r^(ALya)] 

(Note that T = f^'-'^T'.iXi^ya,)/ {l-exp[-T',(Xi^y^)]} for star- 
burst galaxies in the outflow phase.) When q^i <l, continuum 
photons are more significantly absorbed by dust than Lya 
photons, and EWLya can be larger than EW^^^^ (F > 1). In 
our model, although < I, the EW enhancement effect is 
reduced by the dust-independent factor of /q. We plot the EW 
enhancement factor T against the reddening parameter A 1500 
in Fig.[T](right panel) for several values of q^. When q^ < 1, 
r increases with A 1500, and then reaches an asymptotic value 
of fo/qd- In the case of our model (/o = 0.23 and q^ = 0.15), 
the asymptotic value becomes F = 1.53, i.e., only a modest 
enhancement of EW. (It is interesting that a different theoret- 
ical model of Dayal et al. (2008) independently predicted a 
similar value of 1 .6.) Therefore, in our model, we need in- 
trinsically large EW LAEs having EWj^y^ > 160 A to explain 

LAEs having EWLya > 240 A. However, EWlyq is strongly 
suppressed by a factor of F - /„ = 0.23 (or =0.36) in the 
case of A 1500 ^ 0, and hence the dust extinction effect should 
also have an important role to achieve large Lya EWs. 

Our result indicating the importance of dust extinction for 
large EW LAEs is consistent with recent observational results 
that some LAEs indeed seem to have large EWs by the effect 
of clumpy dust distribution (Finkelstein et al. 2008, 2009a, 
2009b). It should be noted that the parameter qj in this work 
is slightly different from an observational estimate of the ge- 
ometry parameter ^obs introduced by Finkelstein et al. (2008). 
The dust-independent factor /o is not taken into account in 

^obs, and it is defined as ^obs = t/'"^ /r^iXLya) assuming screen 
geometry (attenuation cx e~^) in the SED fits. Our definition 
of qa becomes equivalent to ^obs only when we set /o = 1 and 
change the geometry of dust distribution from slab to screen. 

3. COMPARISONS WITH OBSERVATIONS 
3.1. Definition of LAEs and Selection Criteria 

In observations of LAEs, there are mainly two criteria to se- 
lect LAE candidates from photometric samples: the limiting 
magnitude of narrow-band filter that catches redshifted Lya 
lines and the color between narrow- and broad-band filters. 
These roughly correspond to the criteria of Lya luminosity 
(L°b;„ > L%J and Lya EW (EWLy„ > EW* J, respectively. 
Different criteria are applied for different observations, and 
it is important to compare the model with observations under 
appropriate treatments of these criteria. In this work, we al- 
ways select model LAEs by the same threshold values of L^y^ 
and EW^y^ as those adopted in each observation. The thresh- 
old values used in various LAE observations to be compared 
with our model in this paper are compiled in Table [T] 

3.2. Lya Luminosity Functions 

In Fig. |2] we show our model predictions for the Lya LF 
of LAEs, in comparison with observations at z ~3-6. This 
comparison was already done by KTN07 in detail and with 
more observed data, but here we show only the observed data 
that also have UV luminosity and EW data used in this paper. 



As mentioned in the previous section, the thresholds for L*q 
and EWLya are important in this kind of comparison, and here 
we correctly match these conditions between the model and 
the data. Two different panels are shown for the same redshift 
of z = 3.1, corresponding to different observational data by 
different authors using their own LAE selection criteria. On 
the other hand, the LAE selection criteria of the z = 5 .7 data by 
Shimasaku et al. (2006) and Ouchi et al. (2008) are similar, 
and we plot only one model in comparison with them. 

The overall levels and characteristic break luminosities of 
Lya LFs predicted by the outflow+dust (slab) model (the 
model used in this work) are in reasonable agreement with 
the observations. The agreement with the S06 data at z = 5.7 
is rather trivial, because we determined the three LAE model 
parameters by fitting to these data. The bright-end cut-off in 
the model LF for z = 3.1 selected with the criteria of Ouchi 
et al. (2008) is too sharp compared with the data. However, 
according to Ouchi et al. (2008), there is a significant con- 
tamination by AGNs at z = 3.1 and z = 3.7 in the brightest 
luminosity range of L*q ^ 10^^ ergs s"', while no signif- 
icant AGN contamination was found at z = 5.7. This might 
be a possible reason for this discrepancy. Quiescent galaxies 
contribute only to relatively low L*^ LAEs, and their contri- 
bution becomes smaller with increasing redshift. Therefore, 
their contribution is negligible if the limiting magnitude of 
narrow band is shallow (i.e., L*^^ > 10"*^ ergs s"'). 

The alternative models (the simply-proportional model and 
the outflow+dust (screen) model) tested in KTN07 are also 
shown. The difference between the slab and screen dust is 
negligible, while the simply-proportional model overproduces 
the bright-end of Lya LFs especially at lower-z. 

3.3. UV Continuum Luminosity Functions 

Before showing UV continuum LFs of LAEs, we examine 
whether our model correctly reproduces the UV LF of LBGs, 
which are a more abundant and well studied population of 
high redshift galaxies. Because of their simpler selection cri- 
teria based on the Lyman break absorption feature by neutral 
hydrogen in the IGM and ISM (e.g., Yoshii & Peterson 1994; 
Madau 1995), UV LF of LBGs can practically be considered 
as that of all galaxies at the redshift. We show the rest-frame 
UV LFs of all model galaxies (i.e., no selection) at z 3-6 
in Figure [3] Our model are in reasonable agreement with the 
observations. 

We now turn to the UV LFs of LAEs. Figure |4] presents 
comparisons of LAE UV LFs at z < 6 between model pre- 
dictions and observational results. For the LAEs at z = 3.1, 
we plot two different UV LFs for the model LAEs as done 
in Fig. |2] We find that, although the z = 3.1 data points by 
Gronwall et al. (2007) and the z = 3.7 points by Ouchi et 
al. (2008) show quantitative discrepancies (up to a factor of a 
few) at low UV luminosity range (Muv-51og/i < -19 mag) 
in comparison with the model predictions, the overall profiles 
are roughly reproduced. 

3.4. Lya EW Distributions and the Origin of Large EW 

LAEs 

The observed EW distributions by a variety of authors at 
z ^ 3-6 are shown in Figure |5] The statistical la errors and 
upper limits are calculated by the small number Poisson statis- 
tics tabulated by Gehrels (1986). Our model predictions are 
also presented for comparison. A small number of observed 
galaxies have EWLya smaller than the threshold EWl^^,, while 



6 



Kobayashi, Totani, & Nagashima 



the model predicts exactly EWLya > EW*y^. This is because 
the criteria to select LAE candidates from photometric sam- 
ples in actual observations are not exactly the same as those 
in the model based only on Lf^^ and EWf^^. However, the 
difference is not significant and does not affect our main con- 
clusions in this paper 

The model predictions are in good agreement with the ob- 
served distributions at various redshifts. It should be noted 
that our model naturally reproduces the observed EW distri- 
butions well beyond EWLya ^ 240 A, even though our model 
does not include Pop III stars and it assumes the ordinary 
Salpeter IMF. This indicates that Pop III stars or an extremely 
top-heavy IMF are not inevitably required to interpret the ex- 
istence of large EW LAEs. 

To understand the reason why our model could reproduce 
the large EW LAEs, we plot the EW (intrinsic as well as ob- 
servable, as defined in § 12. 5b distributions in the left column of 
Figure|6] for the quiescent, pre-outflow starburst, and outflow 
starburst populations. Here, we show the case of z = 3.1 as 
an example, with no selection about EW (i.e., EW^^^ = A) 
but L*y^ = 10^''^ ergs s~'. The distributions of intrinsic 
EW can be understood by the distributions of metallicity and 
characteristic stellar age of model galaxies, which are shown 
in the middle column of the same figure. Here, the "character- 
istic age" is defined as that of the stellar populations mainly 
contributing to the values of EW, so that EWj^y^^ of the model 
galaxies should be close to those of instantaneous starburst 
having the same metallicity and age. For comparison, the 
contours of EW^^^ of instantaneous starbursts are depicted 
in the same panels, giving an approximate EWj^JJ^ values of 
the model galaxies in this plane. However, the definition of 
the characteristic age is rather complicated, and we explain 
below. 

EW is determined by ionizing luminosity (A < 912 A and 
proportional to the intrinsic Lya luminosity) and continuum 
luminosity around 1216 A, and hence two characteristic time 
scales can be defined: the mean stellar age weighted by ioniz- 
ing luminosity, fj"^, and the mean stellar age weighted by UV 
continuum luminosity, t^^'^. We found that fj^^ becomes sig- 
nificantly greater than fj™ especially at ages of f,™ ^ 10^^ yr, 
and EWj^y^ values of model galaxies are considerably differ- 
ent from those of instantaneous starbursts at either age of fj™ 
or fj^i^. Therefore, we take the harmonic mean"* of the two 
ages to calculate the characteristic age fium, which is the quan- 
tity plotted in the middle column of Fig. |6] Then these plots 
show the effects of age and metallicities in the resultant values 
of EWj"y^ in each model galaxy. 

The intrinsic EW distribution of the quiescently star- 
forming galaxies is narrowly peaked at EWJ"^^ « 70 A. This 
is because these galaxies are forming stars approximately 
constantly on time scales longer than the lifetimes of mas- 
sive stars contributing to the ionizing luminosity or contin- 
uum luminosity around 1216 A, and thus EW reaches an 
equilibrium value determined by stellar metallicity and IMF 
(Chariot & Fall 1993; Schaerer 2003). Although metallic- 
ities of the model galaxies have some scatter in a range of 

We found that the EW of instantaneous starburst population at the har- 
monic mean age is closer to the model EWs, than the arithmetic or geometric 
means. 



Z/Zq ^ 10"' ''-1, equilibrium EWs are narrowly concentrated 
to EWj"y^ « 70 A in this metallicity range. (See Fig. [T]and 
the middle column of Fig. |6l) In spite of the narrow width of 
EW"y^ distribution, EWlyq are distributed more widely be- 
cause model galaxies have a variety of extinction amount, as 
shown in the right column of Fig.|6l Since EWj"y^ is limited 

to be < 100 A, the quiescent galaxy population cannot explain 
galaxies having EWLya ^ 240 A. 

On the other hand, the distribution of intrinsic EW extends 
to a very large value of ^ 400 A in the case of pre-outflow 
starbursts, because of small stellar ages and low metallicities 
(see the middle row of Fig. |6l). The observed EW distribu- 
tion also extends to ^ 400 A, and it is this LAE population 
that contributes to the large EW LAEs of > 240 A. How- 
ever, it should be noted that simply large EWJ^y^ is not suf- 
ficient to have large EWLya. A significant amount of dust 
(i.e., A 1500 > 2 mag) is also required to compensate the EW 
reduction by the dust-independent factor /o. The large EWLya 
found for this population is due to the combined effects of 
stellar population (small age and low metallicity) and clumpy 
dust. 

The bottom panels of Fig.|6]shows the same but for outflow 
starbursts. As shown in the bottom-left panel of Fig. |6l small 
fraction of the galaxies also contributes to the large EW LAEs 
at z = 3 . 1 . They are found to be so young that their EWj^^^ are 
large and that non-negligible amount of dust remains in their 
ISM, which can enhance their EWj^^^ to EWLya > 240 A. 
Hence, the importance of the combined effects of stellar pop- 
ulation and clumpy dust is revealed again. 

3.5. Correlation between Lya EW and Extinction 

The above results suggest that extinction by clumpy dust is 
playing an important role to produce LAEs having large EWs. 
Then we expect some correlations between observable EWs 
and reddening. Therefore we plot EWLya vs. A 1500 (or Ay) 
in the right panels of Fig.|6] The quiescent galaxies populate 
the lower-left corner (small EWs and A 1500), and a clear cor- 
relation between EWLya and A 1500 is found, because EWj^^^ 
has a narrow distribution and the difference of EWLya is due 
to the difference of extinction. The pre-outflow and outflow 
starburst galaxies have a large scatter in this plot, because of 
the wide distribution of EW^^^. However, there is a forbid- 
den region in the lower-right corner, meaning that large EW 
LAEs (EWLya 200 A) must have non-negligible extinction 
(■A 1500 ^ 1 mag). Then we expect some correlation or trend 
between EWLya and A 1500, depending on the relative propor- 
tions of these three populations. 

In fact, recent observations by Finkelstein et al. (2008, 
2009a) have indicated that the extinction by clumpy dust has 
an important effect to produce large EW LAEs at z = 4.5, by 
comparing EWs and A 1200 estimated from the SED fitting. 
The data set of Finkelstein et al. (2009a) can directly be com- 
pared with our model prediction in the EWLya-Ai5oo plane, 
which is given in Fig. [T] (The values of A1200 in Finkelstein 
et al. have been converted to A 1500 by the extinction law as- 
sumed in this work). Direct observational measurements of 
EWLya by Finkelstein et al. (2009a) suffer from large uncer- 
tainties for galaxies with faint UV flux (or, large EWs), and 
hence we also plot model-EWs that are estimated by Finkel- 
stein et al. (2009a) using the UV luminosities calculated by 
the SED fitting rather than observationally measured UV lu- 
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minosities. Because of the degeneracy between dust extinc- 
tion and stellar age, the estimates of A 1500 typically have a 
Icr error of < 0.8 mag (S. Finkelstein, private communica- 
tion). When the three model populations are mixed with the 
LAE selection conditions matched to theirs, we find that our 
model distribution is in good agreement with that of the ob- 
served data. Here, the contribution of the quiescent galaxies is 
negligible because of their relatively large threshold Lya line 
luminosity, L*y^ = lO'*^'^'* ergs s"'. 

The positive correlation among A 1500 and EWLya seems to 
be inconsistent with the results for the LBGs at z ^ 3 (Shapley 
etal. 2003) and those atz ^ 3.5-6 (Pentericci et al. 2009) be- 
cause they reported that more dust extincted LBGs have lower 
EW on average. However, it should be mentioned here that 
their results do not contradict our prediction. This is because 
their samples are limited at relatively luminous range of UV 
continuum, Muv - 5 log/i < -19 mag. In our model, such UV 
luminous galaxies have small EW (EW < 150 A) and less dust 
extinction (A 1500 S 2 mag). The positive correlation predicted 
by our model is not clearly seen in such a small EW and A 1500 
range as presented in Fig. [7] Observations of the galaxies with 
Muv -5 log /z ^ ~19 mag is required to test the validity of our 
model prediction. 

This result gives a further support to the idea of EW en- 
hancement by clumpy dust, which have been inferred from 
the observed EW-A1500 correlation, and independently, from 
the LAE Lya LF modeling in our previous work. Although 
the statistics is still limited, comparisons in this EWlyq-Aisoo 
plane with more data obtained in the future will provide an 
interesting test of our model for LAEs. 

3.6. Mu^j-EW Plot and the Ando Effect 

In Figure|8] the comparisons between the model predictions 
and the observational data in the Muv-EWlyq plane are pre- 
sented. Note that Muv is observable (i.e. dust extinction un- 
corrected) absolute magnitude of UV continuum. The model 
LAEs distribute similarly to the observed LAEs in this plane. 
Particularly, the trend of the Ando effect (smaller EW for 
brighter Muv) is reproduced well by our model. To see the 
reason why our model reproduces the Ando effect correctly, 
we show the EW distributions for several different Muv in- 
tervals in Fig. |9] (Muv is the same as the horizontal axis of 
Fig. [8]) At the brightest Muv, LAEs are predominantly in the 
outflow phase, and their EWs are not larger than ~ 100 A, 
because their dust extinction is not large enough to EW en- 
hancement by the geometrical effect. On the other hand, the 
pre-outflow starbursts that are responsible for the large EW 
LAEs populate mainly the fainter Muv range, because their 
UV luminosities are reduced by extinction. A large amount 
of extinction is necessary to compensate the dust-independent 
factor /o and achieve large EW, and hence a trend of larger 
maximum EW values for smaller UV luminosities appears. 
Therefore, the clumpy dust effect can explain not only the 
existence of large EW LAEs, but also the Ando effect quanti- 
tatively by the same physical process. 

It has been widely discussed that intrinsically larger galax- 
ies have larger extinction on average. For example, Adel- 
berger & Steidel et al. (2000) and Reddy et al. (2008) 
showed that mean extinction increases with the intrinsic (dust- 
corrected) UV luminosity for z ^ 2 - 3 LBGs, and Gawiser et 
al. (2006) indicated that IRAC/Spitzer non-detected (i.e., less 
massive) LAEs at z = 3 . 1 have smaller dust amount than those 
detected by IRAC. This trend may appear to be contradictory 



to our interpretation of the Ando effect, i.e., smaller extinction 
for UV brighter galaxies. However, it should be noted that 
the Ando effect is in terms of the observable (i.e., dust ex- 
tinction uncorrected) absolute magnitude of UV continuum, 
but not in terms of the intrinsic one. There is a considerable 
scatter between the intrinsic and observable UV magnitudes 
because of the strong extinction in the rest-frame UV band. 
As a result, we do not expect a clear trend between extinction 
and observed UV magnitude and hence it is not contradictory 
with the observations. In fact, there are some observations 
that support our interpretation. Reddy et al. (2008) and Buat 
et al. (2009) have reported that the maximum value of atten- 
uation factor by dust increases toward fainter Muv (observ- 
able) while its mean value slightly decreases. In our model, 
the large EW LAEs are produced by strong extinction with 
the clumpy ISM effect, resulting in larger maximum attenu- 
ation factor at fainter Muv- However, the dominant LAEs in 
number density in this magnitude range are, in fact, quiescent 
galaxies with small EW value, as seen in Fig.|9] and hence the 
mean attenuation factor does not increase toward fainter Muv- 

4. DISCUSSION 

4. 1 . LAEs flf z > 6 and IGM Transparency: Implication for 
the Cosmic Reionization 

Although the Lya LF predictions by the KTN07 model are 
in reasonable agreement with the observations at z = 3-6, it 
becomes discrepant rather suddenly with the observed data 
at z > 6. KTN07 argued that this discrepancy can be re- 
solved if the IGM transmission for Lya photons, (the 
fraction of Lya photons that can be transmitted without ab- 
sorption by neutral hydrogen in the IGM), becomes small at 
z > 6. A value of ~ 0.5-0.6 was inferred in order to 
make our prediction consistent with the observed LAE Lya 
LFs at z = 6.56 by Kashikawa et al. (2006) and at z = 6.96 
by lye et al. (2006) and Ota et al. (2008). (Here, is 
defined as a relative transmission compared with that at z < 6, 
and hence = 1 at z < 6. See § 12.21 ) If this interpretation 
is correct, it indicates a significantly higher neutral fraction in 
IGM at z > 6 than lower redshifts, and we may be observing 
the end of the cosmic reionization, although quantitative con- 
version from to IGM neutral fraction xhi = nni/nH is 
rather uncertain (Santos 2004; Haiman & Cen 2005; Dijkstra 
et al. 2007a, b). Here, we examine whether the other statis- 
tical quantities of UV LFs and EW distributions are quantita- 
tively consistent with the above interpretation. 

Figure [TO] shows the comparisons of Lya LF between the 
model and observations at z = 6.56, testing several different 
values of T^^. This is identical to Fig. 4 in KTN07 but here 
for the updated outflow+dust (slab) model. While 0.6 
is preferred in order to fit the bright end of the observational 
Lya LF, model prediction is underestimated by a factor of ^ 2 
at faint end compared with observation. LAE candidates at 
such faint end have not yet been confirmed by spectroscopy, 
leaving room for possible contaminations (N. Kashikawa, pri- 
vate communication). The uncertainty of the observational 
Lya LF at faint end can be as large as a factor of 2 because 
of the low detection completeness of ~ 45% (Kashikawa et al. 
2006). Therefore, in order to avoid the uncertainty caused by 
contaminations and low completeness, we compare our model 
with the other statistical quantities (UV LF, EW distributions, 
and Muv -EWLya plane) by using only photometric samples 
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brighter than L^ya > 10^^ "* h ^ ergs s ' , at which the detection 
completeness is larger than ~ 75%. 

The comparisons are presented in Figure [TT| In our model 
prescription, the Lya luminosities are simply reduced by the 
factor of TJ^™, while Luv is not changed. As a consequence, 
the EW distribution is just shifted into the smaller EW direc- 
tion with decreasing 7]^™, and the positions of model galax- 
ies in the Muv-EW plane simply move downward, while no 
change in UV LPs. However, because of the threshold values 
of L*y^ and EW^^^^, some model galaxies are excluded from 

the LAE selection as T^^^ decreased, resulting in a slight 
change in the faint end of UV LP. This is simply because UV- 
faint galaxies are also Lya-faint on average, and hence they 
have a higher chance to be excluded from the sample by the 
threshold Lya luminosity. 

The EW distribution at EWLya < 120 A shows a better 
agreement with the model when Tl^^ ^ 0.6 than the case of 
^Ly'a' - 1 ■ C)n the other hand, the number fraction of LAEs 
with EWLya 200-260 A seems to favor T^^^ = 1. How- 
ever, it should be noted that the EW measurement is quite 
uncertain for large EW LAEs, because they have faint UV 
luminosity. This can be clearly seen in the Muy-EW planes 
(the right panels of Pig.fTTTl. where we show the vertical dotted 
lines indicating la and 2cr lines of the signal-to-noise of UV 
luminosity measurement. All LAEs with EWLya > 200 A in 
the K06 sample are less than 2a about UV luminosity, and the 
errors of EWLya are quite large, as shown in this plot. There- 
fore, we consider that the EW distribution at EWLya > 200 A 
cannot reliably be compared with the model. When we con- 
centrate on the reliable regions of EWLya < 200 A in the EW 
distribution and Muv - 5 log /K - 1 9 .03 (2f7 line) in Muv-EW 
plane, we find that the model with T^^^ ^ 0.6 gives better fits 

than that with Tj^™ = 1 . Therefore we conclude that the in- 
dication of TI^^ ^ 0.6 obtained by KTN07 is quantitatively 
strengthened by these comparisons. It should also be noted 
that the comparison in UV LP also favors ~ 0.6, by the 
threshold effect about LJ^^'^. 

4.2. Redshift Evolution of LAE UV LF and Selection Effects 

Recently, Ouchi et al. (2008) reported an "anti- 
hierarchical" evolution of LAE UV LPs; they found that the 
LAE UV LPs obtained by the Subaru/XMM-Newton Deep 
Survey show increase of luminosity and/or number densities 
from z = 3.1 to 5.7. Moreover, several observational results 
indicate that the number fraction of LAEs in LBGs increases 
with redshift (e.g., Shimasaku et al. 2006; Ouchi et al. 2008; 
Shioya et al. 2009). This trend has been theoretically quan- 
tified by Samui et al. (2009) as the LAE fraction of ^ 0.1 
and ^ 1 atz = 3-4 and z > 5, respectively, to fit their analytic 
model to the observational data of LAEs. These trends are ap- 
parently in contradiction with our model prediction given in 
Pig. [12 where the UV LP of LAEs continuously decrease to 
higher redshifts beyond z = 3, and the LAE fraction in LBGs 
is approximately constant against redshifts. However, in order 
to discuss such redshift evolutions, it is important to carefully 
consider the dependencies of LAE LPs (both in terms of Lya 
and UV continuum) and LAE fraction in LBGs on the thresh- 
old values of LAE selection: L^y^ and EW*^^. In Pig. [12] we 

have assumed constant values for L}^ and EW^^ regardless 



of redshift. (Note that the LAE fraction in LBGs is not a free 
parameter in our model, but it is derived from more physical 
modelings about fsZ"' taking also into account the adopted 
values of L*y^ and EW*y^ in an observed data set.) Here we 
examine this issue and whether our model is consistent with 
the observations. 

The dependencies of LAE UV LP on the threshold values of 
EW^y^ and Lf^^, predicted by our model at z = 3.1 are shown 
in Pigure [T3] It is found that difference of EW^y^ results 
mainly in changes of the bright end of UV LPs, because the 
UV brightest LAEs have lower EWs on average by the Ando 
effect and hence they are more easily affected by EW^^^. On 
the other hand, L^y^ affects mainly the faint end of UV LPs, 
simply because UV-faint LAEs have small Lya luminosities 
on average although there is a considerable scatter between 
Lya and UV luminosities. In Pigure[T3] the observed UV LPs 
of the LAEs obtained by Gronwall et al. (2007) and Ouchi 
et al. (2008) at the same redshift but with different LAE se- 
lection thresholds are also shown. The behavior of the bright 
end of the observed UV LPs against EWfy^, is quantitatively 
well reproduced by our model. The faint-end behavior against 
Lf^Q, is also consistent with the model prediction, at least qual- 
itatively. 

Our model indicates that the observed anti-hierarchical evo- 
lution of LAE UV LP reported by Ouchi et al. (2008) is 
merely a consequence of a selection effect: decreasing EWfy^ 
with increasing redshift. These results demonstrate the impor- 
tance of the selection conditions of LAEs when the redshift 
evolutions of their statistical quantities are discussed. 

5. CONCLUSIONS 

We have performed a comprehensive comparison between 
a theoretical model of high-z LAEs in the framework of hi- 
erarchical structure formation and a variety of observations 
about statistical quantities such as UV continuum luminosity 
function (UV LP), equivalent width distributions, and corre- 
lation between UV luminosity and EW at redshifts z ^ 3-6. 
The model used in this work was constructed by our previous 
study (KTN07), showing a good agreement with the observed 
Lya LP of LAEs. All the model parameters about LAEs have 
already been determined by fitting to the observed Lya LP at 
z = 5.7 in KTN07, and there is no more free parameters in 
the new comparison made in this work, giving an objective 
test for the validity of our LAE model. In the comparisons 
between the model and observations, we paid a particular at- 
tention about the selection conditions of LAEs; the theoreti- 
cal model predictions are compared separately with different 
observations even at the same redshift, by adopting the same 
conditions as those in each observation for the model predic- 
tions. 

In our model, extinctions by dust for Lya photons and con- 
tinuum photons are treated separately, and the best-fit parame- 
ters to Lya LP obtained by KTN07 indicate that the extinction 
of Lya photons is significantly less effective than that for con- 
tinuum photons around the Lya wavelength. This suggests 
that the dust geometry in ISM of most LAEs is clumpy, and 
this is consistent with recent independent observational stud- 
ies by Pinkelstein et al. (2008, 2009a, 2009b). This clumpy 
dust distribution plays an important role to reproduce the ob- 
servations of the statistical quantities of LAEs considered in 
this paper. 
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We found that the predicted statistical quantities are in 
nice agreement with various observational data in the red- 
shift range of z ~ 3 - 6. In particular, our model naturally 
reproduces (1) the existence of large EW LAEs (> 240 A) 
without invoking Pop III stars or top-heavy IMF, and (2) the 
Ando effect, which is the observed trend of smaller mean Lya 
EW for more UV-luminous LAEs. We carefully examined the 
physical origin of these results, and found that both the stellar 
population and clumpy dust are important to reproduce large 
EW LAEs. EW becomes large when the stellar population is 
young and/or metal-poor, but these effects are not sufficient 
to produce LAEs with EW > 240 A only by these. EW is 
further enhanced by clumpy dust in LAEs with large extinc- 
tion, resulting in large EW LAEs. Our model then predicts 
that there is a correlation or trend that LAEs with larger EWs 
have larger reddening, and we found that this prediction is 
quantitatively consistent with the recent observational results. 
The Ando effect is again explained by the clumpy dust effect, 
because large EW LAEs need a significant amount of extinc- 
tion, and such galaxies have fainter UV luminosities due to 
the extinction. Therefore, we could explain all the statistical 
quantities of LAEs under the standard scenario of hierarchical 
galaxy formation, by normal stellar populations without Pop 
III stars or top heavy IMF, but requiring clumpy dust distribu- 
tion in ISM of LAEs. 

For LAEs at z > 6, the observational data prefer the model 
predictions with a smaller T^^^ than that at z < 6, where 7"^^^ 
is the IGM transmission for Lya photons. The result is con- 
sistent with that of our previous study (KTN07) in terms of 
Lya LFs, but here we confirmed this result by more statis- 
tical quantities (UV LF, EW distributions, and Muv-EWlyq 
plane). Therefore, this result provides a further evidence that 
^Ly'o' rapidly decreasing with redshift beyond z ~ 6, giving 
an important implication for the cosmic reionization. 

The dependence of LAE UV LF and LAE fraction in LBG 
population on the selection criteria of LAEs about Lya lumi- 
nosity and EW was also discussed. We have shown that an 
apparent evolutionary effect can be observed if one uses dif- 
ferent selection criteria for observations at different redshifts. 
It is quite important that the same selection criteria are ap- 
plied to the observational data in order to discuss the redshift 
evolution of LAEs without observational bias. 



There are a lot of theoretical uncertainties in the physics 
of LAEs, and our model may not be unique to explain the 
observed data. However, we emphasize that our work is the 
first to compare a theoretical model of LAEs to almost all of 
available observational quantities (LAE LFs in terms of Lya 
and UV continuum, EW distribution, EW-Muv relation, etc.) 
at various redshifts with a consistent set of model parame- 
ters and with the LAE selection criteria in each observation 
appropriately taken into account. It seems that extinction by 
the clumpy ISM is a unique and necessary ingredient to re- 
produce all of the observed data within the framework our 
model, if we do not consider somewhat exotic assumptions 
such as the contribution from Pop III stars or top-heavy IMF. 
Of course, our model does not exclude the possibilities of Pop 
III or top-heavy IMF, and we need more data and theoretical 
studies to discriminate these possibilities. However, we have 
shown that the standard picture of galaxy formation within the 
framework of hierarchical structure formation can be wholly 
consistent with the available LAE data, simply by introducing 
the clumpy ISM effect for extinction of Lya photons by dust. 

Our theoretical model for the various observational quanti- 
ties of LAEs at various redshifts would be helpful in planning 
an LAE survey at even higher redshifts and interpreting such 
data sets in future studies. The numerical data on these quan- 
tities of LAEs are available upon request to the authors. 
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TABLE 1 

Compilation of Observational Criteria for Lyq EW and Line Luminosity 



z Ref. EW*^[A] logio [L*^/(r2 ergs s-')] /Vpho, Veff [10* Mpc^] 



3.1 G07 20 4 L79 160 3.8 

O08 64 4L69 356 24.0 

3.7 O08 44 42.29 101 23.0 

4.5 D07 14 42.29 ~ 350 50.7 

5.7 S06 20 42.04 89 6.2 

O08 27 42.17 401 31.6 

6.56 T05, K06 7 41.99 58 7.4 



NOTE. — Col. (I): Red.shift. Col. (2): References. T05: Taniguchi et al. (2005), K06: Kashikawa et al. (2006), S06: Shimasaku et al. (2006), D07: Dawson et al. (2007), G07: 



Gionwall et al. (2007), O08: Ouchi et al. (2008). Cols. (3) and (4); Threshold values of rest-frame EWtya and Lya line luminosity, respectively, to select LAEs. Col. (5): Number of 
photometrically selected LAE candidates. Col. (6): Effective survey volume. 
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Fig. 1. — Left: the time evolution of EW|"y'^ ("intrinsic" Lya EW determined by stellar spectra and the case B recombination with the assumption all of LyC 
photons are absorbed and reprocessed into Lya photons in H 11 regions) assuming the S03 population synthesis model with the Salpeter IMF. The solid and 
dashed curves represent the EWj^"^ evolution for instantaneous starburst and constant star formation, respectively. Several curves are shown corresponding to 
different stellar metallicities: Z/Zq = 0, 1 /2000, 1/50, 1 /20, 1/5, and 1, from top to bottom. Right: the EW enhancement factor F as a function of the extinction 
magnitude at 1500 A, A1500, or at V-band, Ay. We plot F for several values of qj (including the value of qj = 0.15 used in our LAE model) in the range of 
qj = 0— 10 as indicated in the figure, while /o is fixed to 0.23 (the value used in the model). 
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Fig. 2. — Cumulative LAE Lyo LPs at z ~ 3-6. The model predictions are represented by the curves, while the observational data are plotted by the filled 
symbols with error-bars. At each panel, the model LAEs are selected with the same threshold value of EWj'j,^ as that adopted in each observation, which is 

indicated in each panel. The vertical dashed line in each panel is the threshold value of L!^^^ of each observation. Top: the predictions from three different /^c" 
model are shown; the simply proportional, the outflow+dust (screen), and the outflow+dust (slab) models are represented by the dotted, dashed, and solid curves, 
respectively. Bottom: the contributions from quiescent, pre-outflow phase starburst, and outflow phase starburst are plotted separately by the dotted, short-dashed, 
and long-dashed curves, respectively, in the outflow+dust (slab) model. The references for data points are Shimasaku et al. (2006), Gronwall et al. (2007), and 
Ouchi et al. (2008). The arrows are taken from Ouchi et al. (2008) which represent the Lya luminosity range where Lya LPs are dominated by LAEs with AGN 
activities. 
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Fig. 3. — Rest-frame UV (A = 1500 A) LFs of all (i.e., including non-LAEs) galaxies at z ^ 3-6. The curves are predictions by our model, while the open 
symbols with error-bars are the observed data of the LBGs at similar redshifts. The contributions from the quiescent galaxies and the starbursts in the pre-outflow, 
outflow, and post-outflow phases are shown separately by the dotted, short-dashed, long-dashed, and dash-dotted curves, respectively. Note that post-outflow 
galaxies are assumed to produce no Lya photon and hence do not contribute to the model LAE population, but contribute to LBGs by UV continuum (see text). 
The references for the data points of LBGs are Steidel et al. (1999), Sawicki & Thompson (2006), Yoshida et al. (2006), Bouwens et al. (2007), and Iwata et al. 
(2007). 
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Fig. 4. — The same as Fig.[3]but for LAE UV LPs. The curves are the model predictions, while the solid symbols with en'or-bars are the observed data of the 
LAEs. The model LAEs have been selected with the same values of L^^^ and EW*^^ with those adopted in each observation shown in each panel (the unit of 
^Lya ^ ~ ^^^^ J^o'h of the top panels are the UV LPs of the LAEs aX z = 3.1, but they are selected with different criteria. The line markings are the same 
as Fig. [3] The references for the data points of LAEs are Shimasaku et al. (2006), Gronwall et al. (2007), and Ouchi et al. (2008). 
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Fig. 5. — Rest-frame EW distributions for LAEs at z ~ 3-6. Tlie curves are predictions by our model, while the histograms are the observations. The observed 
data based on spectroscopically confirmed LAEs are labeled as "spec", while those based on photometric samples as "phot". The error bars and upper Hmits are 
the la Poissonian statistics for small sample number (Gehrels 1986). The line markings of the model curves are the same as Fig. [5] The vertical dashed line 
indicates the maximum EWj^"^ value (240 A) powered by star-formation activity with the Salpeter IMF and the solar metallicity (Chariot & Fall 1993; Schaerer 
2003). The references of the data are given in Fig.|4] except for Dawson et al. (2007). 
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Fig. 6. — The left, middle, and right panels show EW distributions (intrinsic and observable), luminosity-weighted age (iium) versus metallicity, and observable 
EW versus extinction, respectively, of the model LAEs at z = 3.1. The quantity represents a typical age of stellar population that is responsible for Lyo 
EW (see text for exact definition). The contributions from quiescent, pre-outflow phase starburst, and outflow phase starburst are shown separately from top to 
bottom panels. In the left panels, total EW distributions including all of the three populations (quiescent, pre-outflow, and outflow) are also shown by thin curves. 
The small dots in the middle and right panels are model galaxies. In the middle column, the contours of EWj"^^ of instantaneous starburst population (shown in 
Fig.[T) are shown. The model LAEs are selected with = 10^'^ K'^ ergs s"', but with no selection about EWLya. 
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The selection criteria of model LAEs are indicated in the top panel, which are identical to those adopted in Finkelstein et al. (2009a). In the middle and bottom 
panels, the median of the model LAE distribution is shown by the filled circles connected by solid lines, while the 10-90 percentiles are represented by the gray 
regions. The open symbols with error-bars are observational data of 14 LAEs at z ~ 4.5 given by Finkelstein et al. (2009a). The data of A1500 were derived 
by their SED fittings and its typical error is ~ 0.8 mag (S. Finkelstein, private communication). The EWLya data in the middle panel were derived by direct 
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Fig. 8. — Distributions of LAEs in tlie Muv-EWLya plane at z ~3-6. In tlie upper panels, model LAEs are plotted by small dots. In the lower panels, the open 
symbols are corresponding to each of the observed galaxies. The filled circles connected by solid lines are the mean of the model galaxies, and the gray region 



represents the 10-90 percentiles of the model galaxies. See Fig.|4]for the references. The solid curve and horizontal line in each panel indicate Lf]^^ and EW| 



Lya ' 



respectively. Note that M\]\ is in the rest-frame 1500 A, and we must convert it in the rest-frame Lya wavelength to depict the curve of Lf^^^- We assumed a SED 

with /3 = — 2 as an approximation here. The dotted curves show the contours of Lya line luminosity corresponding to 2 X 10*^^, 10*^^, and 5 X 10*~ ergs s" 
from top left to bottom right, assuming /3 = —2. 
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Fig. 10. — Comparisons between the cumulative Lya LF data of Kashikawa et al. (2006) at z = 6.56 and model predictions with three different values of the 
IGM transmission for Ly« photons, T^^^ =10 (top), 0.8 (middle), and 0.6 (bottom). The thick solid curve is for all model LAEs, and the other curves are for 
model LAEs with different phases (the same line markings as those of the lower panels of Fig.|2). The filled squares are all LAE candidates in the K06 sample, 
while the filled circles are those with spectroscopic confirmation. The vertical dashed line indicate s th e threshold Lya luminosity of lO''^'* If^ ergs s"', which 
will be applied in the comparisons with UV LF, EW distribution, and Myy -EWLya plane (see Fig. II It . 
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Fig. 11. — The same as Fig. 1101 but for comparisons with UV LF, EW distribution, and Muv~EWLya plane from the left to the right panels. The line 
markings of the model curves in the left, middle, and right panels are the same as those of Figs.|4]|5] and[8] respectively. Here, only LAE candidates having 
I^"^ > 10''^ * If^ ergs s"' are used from the K06 sample, and this condition is included in the model predictions as well. The vertical dotted lines in the UV LF 
and Muv ~EWLya plane indicate the 1 and 2 a limiting UV continuum magnitudes of Kashikawa et al. (2006). 
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Fig. 12. — Redshift evolution of LAE UV LF (left) and LAE fraction in LBGs as a function of UV luminosities (right) predicted by our model from z = 3.1 
to z = 8.8 with = 1. The dashed curves in the left panel show the UV LFs of all galaxies (or observationally, LBGs). The galaxies with L\ 
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h ergs s and EWLya > EW* = 20 A are selected as LAEs regardless of redshift. 
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Fig. 13. — UV LFs of LAEs at z = 3.1 predicted by our model with different values of threshold EW* ^ (left) and L'^^^ (.right) in LAE selection. The thin 
solid curve represents the UV LF of all galaxies including non-LAEs. The observed data of UV LFs of LAEs at z = 3. 1 obtained by Gronwall et al. (2007, G07) 
and Ouchi et al. (2008, O08) are also shown. Their LAE selection thresholds are ll%f„/(h-^ ergs s"'), EWj^y^/A] = (Iff*' ™, 20) and (Iff*'-*', 64) for G07 and 
O08, respectively. 



